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ABSTRACT 

We study the classical motion in bidimensional polygonal billiards on 
the sphere. In particular we investigate the dynamics in tiling and generic 
rational and irrational equilateral triangles. Unlike the plane or the neg- 
ative curvature cases we obtain a complex but regular dynamics. 



1 INTRODUCTION 



In this work we consider the classical motion of a point particle inside a two-dimensional 
polygonal billiard on a surface with constant positive curvature. Flat billiards and 
billiards on a surface with negative curvature have been extensively studied. It is well 
known that the dynamics in flat polygonal billiards depends on whether the inner 
angles of the polygon are rational multiples of n |], |], |3|. If this is the case these 
systems are referred to as 'pseudointegrable' since they possess two constants of mo- 
tion and the flow is restricted to a two-dimensional invariant surface. If at least one 
of the vertex angles is an irrational multiple of 7r the polygon is believed to possess 
ergodicity @]. On the other hand the interest in studying polygonal billiards on a 
surface with negative curvature is that the classical motion is simple and as chaotic 
as possible since the flow on these surfaces is hyperbolic 0. 

By investigating the dynamics on a spherical surface we explore the other limit: in- 
stead of having more chaoticity, more focusing is expected with respect to the planar 
case as a consequence of the positive curvature. The motivation for studying these 
systems is to see how this focusing mechanism together with the compactness of the 
surface affects the classical motion. In particular we will investigate the condition for 
integrability and stress the importance of tiling billiards. In the general case we will 
explore numerically the phase space portrait of curved polygonal billiards (i.e., peri- 
odic orbits, invariant surfaces, etc), and try to explain these numerical observations 
in a rather intuitive way. 

The paper is organized as follows. In section 2 we present the formalism. In section 
3 we briefly discuss the case of billiards enclosed by meridians and parallels and then 
concentrate on the study of polygonal billiards, i.e. billiards whose boundary consists 
entirely of arcs of geodesies. In particular we will focus on equilateral triangles. 
We will first study tiling triangles in subsection 3.1 and then generic triangles in 
subsection 3.2. For this we use the entry-exit map and finally make a comparison 
of small curved triangles with triangles on the plane. Conclusions are presented in 
section 4. 

2 THE MODEL 

The motion on the sphere will be described in terms of the spherical coordinates: 
the polar angle a and the azimuthal angle f3. The coordinate curves a = const and 
(3 = const form an orthogonal net. The line element for a sphere of radius R has the 
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usual form: 



ds 2 = R 2 da 2 + R 2 sin 2 adj3 2 (1) 

and the curvature is 1/R 2 . 

The geodesies are the great circles of the sphere and can be viewed as its intersection 
with a plane passing through the origin. They are labeled by the coordinates (9 G , <p G ) 
of the unit vector normal to this plane and their equation reads: 

cot# G 

tana = — (2) 

cos(,9 - <j) G ) 

Note that a and (3 are coordinates on the physical sphere where the billiard lies while 
9 G and 4> G denote a point on what we will call the dual sphere. This is illustrated in 
Fig.l. In the problem we are considering the particle moves freely along an oriented 
geodesic labeled by (9 G , (f> G ) until it suffers a specular reflection at the boundary and 
jumps on a different geodesic. Point (9 G , (p G )i maps into a new point (9' G , 4>' G )i+i 
and a given trajectory can be specified uniquely by an infinite succession of pairs 
(...(9 G ,<f) G )j...) on the dual sphere labeling the successive geodesies on which the 
motion takes place. In this way the motion is completely described by a map acting 
in discrete time on a bounded domain on the surface of the dual sphere. The area 
element dA = sin 9d9d(p is preserved by this bounce map T. Thus cos9 G and <p G 
are canonical coordinates on the dual sphere. They provide a convenient reduced 
description of the classical motion alternative to the Birkhoff coordinates. 



3 POLYGONAL BILLIARDS 

The simplest systems on the sphere one can think of are convex billiards whose 
boundary consists of meridians (f3 = const) and parallels (a = const). These are not 
polygons since parallels (except for the equator) are not geodesies. It is easy to see that 
in reflections on meridians and parallels, the quantity sin 6? is conserved . Therefore, 
the dynamics for these billiards is integrable. Depending on the initial conditions 
we get periodic and non periodic orbits, giving rise to rational and irrational tori in 
phase space. 

We now concentrate on polygonal billiards, i.e. systems enclosed by arcs of geodesies. 
Note that in this case not only the trajectories but also the boundaries can be specified 
by points on the dual sphere. Therefore the problem can be viewed as a mapping of 
unit vectors by a matrix also expressed in terms of coordinates on the dual sphere. 
This would not be possible in the general case of a billiard with arbitrary enclosure. 
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For simplicity we will restrict our discussion to the case of equilateral triangles. While 
flat equilateral triangles are integrable, spherical equilateral triangles present a rich 
variety of possibilities depending on their size. This is due to the existence on the 
sphere of a definite relation between the size and the shape of the triangle. The inner 
angle u can take any value in the interval [|,7r], and the area (that can be related to 
the total curvature) is: A = R 2 (3u — n). 

For a triangle centered on the north pole the three sides will be specified by the 



vectors 



lB ±B\ 



which label the intersecting geodesies. Here, 



1) 



f 1 with 



1, 2, 3, and the angle 9 B is related to the inner angle uj through: 



sin0 J 



cos I 

COS I 



(3) 



Vertex Vi defined as the intersection of geodesies (i — 1) and i will be located at 
(a v ,$Y) , where (a v , (3^) denote coordinates on the physical sphere given by: 



tana 



v 



tan 9 B 



(4) 



and 



=*f + ^ (5) 

It will be useful in the following to introduce for each vertex Vi the curve Cj, locus of 
the points (#,</>) specifying all the geodesies passing by Vi. Its equation reads: 



2tt 

+ t 



v 



tan^,' = — - 



cot a 

cos(0, - (3Y) [ J 

It follows from Eq.(||) that curve Cj is itself a geodesic on the dual sphere labeled by 
(a v , PY). Note that this is not true in general but rather a consequence of the metric 
of our problem that will play an important role when understanding the structure of 
phase space. Once we have defined the discontinuity curves Cj it is easy to see that a 
point (9, 0) corresponds to a geodesic entering (exiting) side i if: #;(0) < 9 < 9 i+ i((p) 
( 9i+i(<p) < 9 < 9i((f>) ) . 

The area preserving mapping T can be written explicitely as an orthogonal matrix 

/ sin#cos</> \ 



acting on the unit vector 



sin 9 sin < 

cos 9 



as: 



/ 



T 



cos< 



sm 



sm i 



\ / -cos26» B sin 29 B \ ( cos 



lB 



kB 



COS 



kB 



V o 





1 / 





sin 29 B 



-1 
cos2# B J 



— sini 




sm i 

cos 



(7) 
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for 9i(4>) < 9 < 8i + i((f>), that is for a geodesic entering side i. 

Since the problem is symmetric under rotations in ^ around the center of the triangles 
it will be useful to define the operators T + and T~ as : 



T 



± 



/ cos f ± sin ^ \ / - cos 26 B sin 26 B \ 
Tsin^ cos^f ■ -1 



x 3 I J 



(8) 



V sin2# B cos 20 s / 

With this definition according to which, for example, T 3 T 2 Ti = T + T + T + and T 2 Ti = 
T + T~, we indicate if the particle exiting from a given side hits the next or the previous 
side in increasing order with i. 

The trajectories in the triangle will be classified according to an infinite symbol se- 
quence obtained by listing the sides successively hit. The code alphabet will consist 
of +, — signs. This classification, as we will see, is not one to one in the sense that 
even infinite length sequences do not distinguish single trajectories uniquely. 

3.1 TILING TRIANGLES 

We now study the case of equilateral triangles that tile the sphere under the reflection 
rule. Since the rotation group has a finite number of discrete subgroups there are 
only a few ways of tessellating the sphere. Tiling triangles are such that their vertices 
coincide with those of a face of a regular polyhedron. The three possible cases are 
the tetrahedron with ui = -y, the octahedron with u — | and the icosahedron with 
u = ^r-. We will see that tiling triangles constitute a very particular class of billiards 
that not only are integrable but for which only periodic orbits are present. 
To study the motion in these three particular triangles we follow the procedure pre- 
sented in HI : every time an edge is hit, instead of reflecting the incident geodesic, 
we reflect the billiard across the edge and follow the same geodesic into the replica. 
That is, the motion is viewed as a unique geodesic entering and exiting copies of 
the original billiard. Since the surface is compact and the billiard is tiling it is clear 
that only periodic orbits exist. In order to determine their periodicity n p we have 
to fold back, for each copy, the corresponding segment of geodesic into the original 
billiard. When, after a certain number of these operations, the image coincides with 
the original the trajectory closes. To do this we have to consider the symmetry group 
of the corresponding polyhedron and the possible circuits on it, that is, which faces 
are visited by the geodesic. 

Triangle with u> = ^ . 
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Two circuits are possible on the projected tetrahedron: one visits 3 faces, correspond- 
ing to orbits of type [... +—...] in the triangular billiard, and the other 4, corresponding 
to orbits [... + + + ...] . Both circuits are indicated in Fig. 2, where the notation for 
the tetrahedron is defined. Sides are indicated with a, 6, c..., while with cr a , crj,, u c ... 
we denote reflections across these sides. Face / corresponds to the original billiard. 

The orbits of the first type correspond to following a geodesic passing through 
faces : /, a^I, <t c <T{,/, cr a i7 c abl, .... The transformation <r j&i is a rotation Cf through 
about the line of intersection of the reflection planes. Thus the product of 6 
reflections, corresponding to a rotation in 47r, brings back the segment of geodesic 
to its original position and orientation. Therefore, general periodic orbits of type 
[... H — ...] have periodicity n p = 6. 

There are three particular orbits of lower periodicity n p = 2. These are the ones 
invariant under the rotation around three of the four C3 axes of the tethahedron, 
that is, the geodesies resulting from the intersection with the sphere of the planes 
orthogonal to these axes. 

The orbits of type [... + + + ...] correspond to a geodesic passing through faces: 
I,<7bI,o'cO'bI , OjU c Ubl ,a a afa c abl , ■■■■ It can be seen that the product of three re- 
flections OkOjOi is a rotary-reflection S4 about a one of the binary axis joining the 
centers of two opposite sides of the tetrahedron, as shown in Fig. 2. Since (S4) 4 = E, 
12 reflections are needed in order to recover the initial segment of geodesic in the 
original billiard and the periodicity of the orbits will be n p = 12. There is, finally, an 
orbit with periodicity n p = 3 that joins the centers of the tree sides of the triangle. 

Triangle with uj = | . 

Due to the geometry of the octahedron, all the trajectories in this triangle are 
of the type [... + + + ...]. A product of three reflections is equivalent to an inversion, 
transforming 9 G , in n — 9, tp + n. Since any geodesic is invariant under this operation, 
the periodicity for all the orbits will be n p = 3. 

Triangle with uj = ^ . 

Following the same procedure we find in this case orbits of type [... + + + ...] 
which have in general periodicity n p = 15 (and one with n p = 3 that joins the centers 

of the three sides of the triangle) and orbits of the type [... + H ...] with n p = 12 

(and three orbits with lower periodicity n p = 4). 

One can also think of triangles that are not tiling in the strict sense but, by successive 
reflections on their sides, cover the sphere more than once. This is the case of a triangle 
with uj = that covers the spherical surface twice. It has periodic orbits of type 



5 



[... H — ...] with periodicity n p = 10, of type [... + + + ...] with periodicity n p = 15 and 
of type [... + H h H h H — ...] with periodicity n p = 9. 

3.2 GENERIC TRIANGLES 

Let us now investigate the case of a generic triangle with an arbitrary 9 B . Following 
@ we will divide the available phase space in the (9, 0)- plane in entry (or exit) 
domains of the three sides. The enter (exit) domain is defined as the set of points 
associated with the oriented geodesies that enter (or exit) a given side. In order to 
do this partition, we use curves C iy defined in Eq.flU). 

In FIG. 3a we show an example of an entry-exit map. The intersections of curves 
Ci and Ci-i correspond to the points associated to side i — 1 with the two possible 
orientations. Each point in the available phase space belongs to one entry and one 
exit domain. Each exit domain is intersected by two entry domains. 
After n iterations of the map according to Eq.(^) the phase space gets partitioned 
into domains which are enclosed by geodesies, since the images by reflection of the 
separating curves Cj are also geodesies. Each domain can be labeled by a sequence 
of n symbols denoting the n sides successively hit by the trajectory. 
The only general rule limiting the possible sequences is that a side cannot be hit twice 
consecutively. But it is clear that for each particular triangle the existence or absence 
of a region associated with a given sequence is determined by the geometry, in this 
case, the inner angle u. For example the domain corresponding to an infinite periodic 
sequence of type [... H — ...] exists only in triangles with angle uj > |. 
FIG. 3b shows the partition of phase space after two iterations of the map. 
As long as the number of iterations remains finite the different allowed domains are 
bounded by arcs of geodesies and therefore polygons of increasing complexity. When 
a given sequence of symbols is repeated periodically an infinite number of times the 
resulting set of trajectories is a chain of elliptic islands bounded by a smooth curve 
which is the limit of the above mentioned polygons. Short repeated sequences lead 
to large and fairly regular areas, while long sequences lead to complex island chains. 
In Fig. 3c shows the domains corresponding to the repetition of some short codes. 
Before studying the structure of these islands in more detail, let us have a look at 
the entry-exit maps for the special case of the tiling triangles treated in the previous 
section. The common feature to these particular triangles is that the corresponding 
separating geodesies Cj are reflected into themselves after a few iterations. Therefore 
phase space gets partitionned into a few polygonal (in the sense that they are enclosed 
by a finite number of segments of geodesies) domains corresponding to different peri- 
odic codes. Inside each domain, all orbits are periodic with the same periodicity and 
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we recover the result obtained in the previous section. This is illustrated in FIG. 4. 
Other rational but non-tiling triangles may have domains in which all trajectories are 
periodic, coexisting with the families of elliptic islands. These triangles with periodic 
domains can be found by requiring the angle 9 B to be such that successive appli- 
cations of the bouncing matrix corresponding to a chosen sequence gives the unit 
matrix. For example requiring that (T + T + T + ) n = I gives 

qos — 

9 B = arccos £- (9) 

cos I v ' 

o 

while (T + T~) n = I fixes 

cos 2 

# B = arcsin (10) 

cos I K ' 

where k and n are integers. 

Thus, if a triangle defined by 9 B and its complement, i.e., the one defined by | — 9 B 
have inner angles uj\ and uj2 which are both rational, angle 9 B satisfies simultaneously 
Eq.(H) and Eq. fllOl) and the pair of rational angles (ui,^) fulfill: 

coscui + cosa;2 = — - (11) 

We found numerically two pairs of rational angles satisfying Eq.(|ll]): {u)\ = f,^2 = 
and (ui = ^f,0J2 = 4f), corresponding to the triangles analyzed in the previous 
section. 

Summarizing we see that the entry-exit map is a useful tool to determine whether, 
according to the geometry (in this simple case the inner angle u), a domain corre- 
sponding to a given code exists or not in a given triangle. For this it is sufficient 
to reflect the segments of the separating geodesies C% enclosing the initial domain 
according to the chosen sequence. The image domain after infinite iterations can be 
either empty or finite, in contrast with the hyperbolic case where non empty domains 
associated to a code are reduced to a single point. Finite domains are in general 
bounded by smooth curves, resulting from the intersection of an infinite number of 
geodesies. In rational triangles some particular codes are associated to polygonal 
domains enclosed by a finite number of geodesies. 

In order to understand the detailed structure of phase space in the generic case we 
now follow individual trajectories. A phase space plot is shown in Fig.5 for different 
initial conditions. It reminds very much the stable regime of a sawtooth map [0. 
At the center of each island, belonging to a set of n p islands, sits a periodic orbit of 
period n p surrounded by a family of nested invariant curves. These correspond to 
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open orbits having an infinite but periodic sequence of period n p . That is, a given 
code does not specify a single trajectory, as in the hyperbolic case, but an entire 
family of orbits. The size of the islands decreases as n p increases and phase space 
takes a fractal structure. 

The existence of these islands is a consequence of the focusing mechanism on the 
spherical surface. The fact that two geodesies on the sphere intersecting at Ointi&int 
cross again at n — 6 int , (p int + 7r reflects on the structure of the orbits. To illustrate this 
we consider a periodic orbit with initial conditions (s pM ', P\\°) expressed in Birkhoff 
coordinates and an open trajectory close to it. If the initial conditions are close enough 
to the ones of the periodic orbit the particle on the open trajectory will follow the 
sequence of the periodic orbit and, after a period, end at a point s on the boundary 
that is shifted with respect to the point s p -°\ The position s as a function of the 
number of period iterations is plotted in Fig. 6 for a given orbit in a curved and a 
flat triangle of equal area . In the flat billiard the shift s — s p o ' increases linearly, 
in such a way that, after a certain number of iterations, the vertex will eventually 
be reached and the original sequence broken. In the spherical case the shift is an 
oscillatory function of the number of iterations. If the amplitude of oscillation is 
small enough the corner is never reached: the pendulating motion lasts forever and 
the particle keeps repeating the original symbolic sequence but never retracing itself 
exactly. This is what happens in the islands of Fig.5. As long as the concentric curves 
do not become tangent to curves Cj, that act as separatrices, the open trajectories 
have the same code as the central periodic orbit. When the outer curve gets tangent 
to Ci, the original sequence is broken and a new and more complex sequence appears, 
corresponding to a new chain of islands. 

The corresponding transformation is a product of a string of orthogonal matrices Tj. 
The resulting 3x3 matrix is itself orthogonal and therefore has an eigenvalue 1, cor- 
responding to preservation of the norm, and a pair of complex conjugate phases e ±tn . 
Thus the motion inside each chain of islands is labeled by a code and is everywhere 
elliptic with the same rotation angle Q. This is clearly displayed ing Fig.5 and we 
have verified it numerically for several short sequences. 

An alternative way to analyze the phase space structure of triangular billiards on the 
sphere is to go to the limit of small triangles where the total curvature can be seen 
as a perturbation from the integrable case of flat equilateral triangles. This limit 
corresponds to take u — ► |, (that is 9 B = | — e). This restricts the available space 
phase to geodesies with 9 G close to f . For the segments inside the billiard, that is (3 
in a small interval around G + n, the equation of the geodesies reads: 
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where = | — 6 = e and $ — <f> + it, that is, the equation of a straight line in 
polar coordinates in the plane. 

The reflection across a side of the planar triangle is given by: 

a' = - a + 26 cos(p - <f) (13) 

P' = -P + 20o (14) 

which can also be obtained by expanding Eq.(|7|) to first order in e. All trajectories 
lie on tori filling phase space. In Fig. 7a several rational and irrational tori are shown. 
Under perturbation, these rational tori are destroyed. In contrast with the generic 
case for which elliptic and hyperbolic cycles of fixed points survive, fixed points that 
survive perturbation in the spherical case are all stable. The absence of hyperbolic 
trajectories might be attributed to the fact that Lyapunov exponents are imaginary 
as a consequence of the focusing mechanism due to the positive curvature. This point 
needs further investigation. 

Around each stable point of the periodic orbit n p a family of invariant curves, cor- 
responding to open trajectories, develops, (that is, we get a chain of n p islands). 
This is shown in Fig. 7b. As mentioned in the previous section, the periodic code 
corresponding to the central periodic orbit characterizes the whole family. 

4 CONCLUSIONS 

This work was a first step in the understanding of the classical dynamics in billiards 
on a spherical surface. We only considered polygonal billiards and, more specifically, 
equilateral rational and irrational triangles. The main observation was the existence 
in the vicinity of each periodic orbit of a set of open orbits characterized by the 
same infinite periodic code. This is a consequence of the focusing mechanism and 
constitutes a substantial difference with respect to the planar and hyperbolic cases 
where each infinite code is associated to a single orbit. The map has a regular but 
very complex structure. The phase space is covered by chains of elliptic islands, whose 
multiplicity increases and size decreases with the period. Tiling billiards constitute 
a very special system where only periodic orbits are present. We have also studied 
non-equilateral triangles and the general structure of phase space is very similar. We 
do expect differences if the boundaries are not geodesies: in particular the motion 
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is likely to develop large chaotic regions. This question together with the issues of 
quantization will be addressed in the future. 
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(a) (b) 

Fig.l - (a) shows three intersecting geodesies enclosing a triangular billiard centered 
at the north pole. In (b) the shadowed area is the corresponding phase space on the 
dual space to the geodesies. 




Fig. 2 - (a) shows schematically a circuit visiting 3 faces of the tetrahedron, corre- 
sponding to an orbit of type [...+ -...] in triangle I. In (b) a circuit visiting 4 faces, 
associated to an orbit of type [...+ + +...]. 
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Fig. 3 - (a) Entry-exit map for a triangle with 9 B = 0.7. (b) Phase space after two 
iterations gets partitionned in domains corresponding to orbits of type [+ +] ([- -]) 
and of type [+ -] ([- +]). (c) shows the domains corresponding to orbits [...+ -...] 
(1), [...+ + +...] (2) and [... + + + H ...] (3) after infinite iterations. 
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Fig. 4 - Phase space portrait for tiling triangles, (a) For u = | curves Cj reflect 
into themselves, (b) For to = ^ curves Cj reflect into C {6 B = |). (c) For a; = 
curves Cj reflect into C[ which, in turn, reflect again into Cj. 
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Fig. 6 - Position on the boundary s plotted as a function of the number of iterations 
for an open trajectory in the neighborhood of a periodic orbit of period n p = 3, 
starting at s pM ,. Circles correspond to the curved triangle, diamonds to the planar 
triangle of equal area. Horizontal lines indicate the edges. 
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Fig. 7 - (a) Some tori in a planar triangle. Circles correspond to a rational torus 
with n, p = 3, squares to n p = 34 and the straight line to an irrational torus, (b) 
For a curved triangle with cu = ^jjf- open trajectories develope around the elliptic 
points. 
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